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Equations to investigate fluctuations in cosmic X-ray background radiation due to 
point-like sources at high-redshift are formulated in a systematic way. The angular power 
spectrum of X-ray background fluctuations is investigated from large-scales to small- 
scales in various cosmological models such as open universe models and models with 
the cosmological constant, assuming a simple evolution model of the sources. The effect 
of epoch-dependent bias is demonstrated for small-angle fluctuations. The contribution 
from shot noise fluctuations is also discussed. 
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. I. INTRODUCTION 

4^ ■ 

Cosmic X-ray background (CXB) has been studied for long time in astrophysics |5J]. Recently several 
authors have discussed the possibility of the CXB fluctuations as a probe of the structure formation 
in the high redshift universe J^-§|. These works are motivated from observational results of the deep 
survey carried out by ROSAT, which resolved a significant fraction of the X-ray background into sources 
(~ 70%). It turns out to be that most of the sources are extra-galactic objects, i.e., AGN and other 
X-ray luminous galaxies in the high redshift universe. While situation is not so clear in the harder band, 
something similar would be happening. Several X-ray mission projects are in progress H, which will give 
us a fine solution of the source problem in near future. 
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On the other hand, a scenario of the cosmic structure formation is now constrained from the observa- 
tions, e.g., temperature anisotropies in the cosmic microwave background radiation (CMB) and large-scale 
distribution of galaxies. Future satellite experiments of CMB and survey projects of the large-scale dis- 
tribution of galaxies will provide severer constraints on theoretical models of the large scale structure 
formation. At present the scenario introducing Cold Dark Matter (CDM) seems to be successful though 
some modification may be required |l0| , [ll|] . 

Formation processes of subgalactic objects, galaxies, AGNs and clusters at the high redshift universe 
are one of the hottest topics in the area of astrophysics and cosmology. Since these objects are parts of 
the large scale structure of the universe, we could expect that information about the structure formation 
in the high redshift universe would be obtained by investigating them. 

Several works have been done to answer following questions: Which kind of information about cosmic 
structure formation can be obtained from an analysis of the X-ray background fluctuations? Are they 
really a probe of cosmology? For example, Boughn, et al. claimed that the cross correlation between 
the X-ray background fluctuations and the microwave background anisotropies can provide a constraint 
on a cosmological constant. Small-scale fluctuations in the X-ray background have been investigated in 
njl2-|l4|. On the other hand, Lahav et al. calculated the angular power spectrum of the fluctuations 



by a multipole expansion method. They mainly focused on large scale fluctuations and their formula is 
only limited to the case of the flat universe. Recently Treyer, et al. have compared the theoretical models 
by Lahav et al. with the observational data of HEAOl ||. 

In this paper we expand Lahav et al.'s formula to more general cosmological models, i.e., open universe 
models and models with the cosmological constant, and investigate general behaviors of the angular 
power spectrum of fluctuations from large-scales to small-scales. The X-ray background fluctuation in 
the open universe has been considered in ref. |I5| ], and the similar expression was obtained. However, in 
that paper, the discussion was only focused on the large-scale fluctuations due to the source clustering, 
and the derivation of equations was very complicated. In this paper we formulate the treatment in a 
simple way. And the shot noise fluctuation is treated in a systematic way. The formulation developed in 
this paper will be useful for studying the clustering of high redshift sources not only for X-ray sources 
but also for other point-like sources such as radio galaxies [[l6 

This paper is organized as follows. In section 2 we derive expression for the angular two-point correlation 
function by using the multipole expansion method in a systematic way. In section 3 the fluctuations in 
the X-ray background are investigated assuming a simple luminosity function of X-ray sources. The 
competition of the fluctuations due to the source clustering and the shot noise is also studied there. 
Section 4 is devoted to summary and discussions. Throughout this paper we use the unit c = 1. 



17 and 7-ray bursts. 



II. FORMULATION 



In this section we formulate equations to calculate background fluctuations which come from point-like 
sources at high redshift. Essence of the formulation is divided into two parts. One is how the sources 
are distributed in the universe. This is the problem of source evolution and statistics of the distribution, 
which we describe in the latter part of this section. The other part is how the sources are observed in 
the sky. Which is obtained by solving light propagation in the expanding universe, if the distribution of 
the sources in the universe was given. 

First we consider the light propagation. The photon intensity I v in an expanding universe is described 
by the radiative transfer equation: 



drj dxi 



dv 



dl v 
drj d^i 



77), 



(2.1) 



where a is the scale factor which is normalized to be unity at present, 77 is conformal time defined by 
adr\ = dt with the cosmic time t, H. is the Hubble parameter defined as a/ a with over-dot denoting 77 
differentiation, 7* is the directional vector of the photon momentum, and j v (x) is a field of emissivity per 
unit comoving volume. Assuming that an (i)-th point-like source is located at the coordinate and 
has the power put out per unit frequency and per unit time Lu {if) , we write the field of emissivity as 



j„(x,„) = — ^TlM^^x-xW). 



(2.2) 



We solve equation (2.1) with (2^2) in the Friedmann- Robertson- Walker space-time with the line element 

ds 2 = -dt 2 + a\ j dx i dx j , (2.3) 
where 7^ is the three-metric on a space of constant negative curvature K: 
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(2.4) 



In order to take the limit of the flat universe, let us first introduce a radial coordinate r instead of \ by 
X = V ' —Kr and take the limit K — > 0. The scale factor a is governed by the Friedmann equation, 
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H 



h$[ -n + n K + a 2 n A ) , 



(2.5) 



where Hq is the Hubble parameter with Hq = 100 h km/s/Mpc, f2o is the density parameter, J7a = A/3Hq 
which is the density parameter of the cosmological constant A, and Qk = 1 — — which describes 
the spatial curvature of the universe. 



Employing a new variable / = I v jv\ we can rewrite equation (2.1) as 



df | • df Uv dJ_ | d 7 J df _ a 
drj dxi dv drj drf 4Tr(av) 3 



£4^)<5( 3 )(x-x«). 



Introducing vq{= va) instead of v 1 (2.6) reduces to 
df , 7 9/ , drf df 



£^/.fo)* (8) (*-* 



Operating J"^ 2 dv^VQ 3 on the both side of equation (2.7), we have 
dl 



drj 



' dxi drj dji 4tt 4^ .L "/-W^ 



where J = J^ 2 di/QVo 3 f ' . Since equation (2^ 



is rewritten as 



d?7 3?7 d?7 <9 7 J 47T-4-' v ; wy w/ 



where the luminosity of the i-th source is defined by 



and K^(r/) is defined by 



(0f«1 - 



dvc\L 1 (77), 



then by integrating (2.9), we get 

/(t7o,x ,7) = ^E / ^(l)^ (i) (l)i (i) (# (3) (x(l,7)-x (i) ), 



(2.6) 



(2.7) 



(2.8) 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



where x(77, 7 ) stands for a photon path. 

Because £\ #( 3 ) (x — xW ) is regarded as a number density field and we replace it with 



^V 3) (x-x w ) => J d\ogLn(L,r],x) 



(2.13) 



* K 1 (rf) / a(rf) denotes the usual K-correction. 
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where n(L, 77, x) is the luminosity fun ction which denotes a number of sources per unit comoving volume 



and per unit logL. Then we rewrite (2.12) as 



7(r7o, x , 7) = J dlogLL J dr]a(r])K(L,ri)n(L,r},x(T],^)), (2.14) 



where K(L,r\) is defined for a source with luminosity L in s ame w ay as (2.11). 



Here we mention the range of 77 integration in equation (2.14). The range of 77 integration depends 
on observational situation and strategy. In this paper we focus on background fluctuations. In order to 
reduce the shot noise fluctuations, which are described later in this paper, we assume that bright nearby 
sources are removed from an observed map. Darker and darker sources we get rid of, more and more 
distant sources are removed. This flux cutoff limit is denoted by S c . The local luminosity of a source L 
is related with the observed flux S and the distance, which we describe by the conformal time 77 in stead 
of the redshift, as 

aLK(L, 77) 

S = toDW ' (2 ' 15) 

where -D(?7) = (— K)^ 1 / 2 sinh(\/— ^(770 — v)) an d Vo is the conformal time at the present epoch. Note 
that the distance D is related with the luminosity distance by D(rj) — + z) flS]. In case of 

observations with the flux cutoff limit S c , the range of 77 integration must be < 77 < 77 c , where 77 c is 
determined by solving ( 2.15| ) with setting S — S c . Thus ?7 C is a function of L and S c , in general. As 



we will see in the next section, the removal of bright sources is an important factor in predicting the 
fluctuations. 

In order to discuss statistics of the fluctuations we define the angular two-point correlation function: 
°( )= J J ^ ^J 2 *(7i ■ 72 ~ cos 0)/(r7o, x , 71)^0, x , 72) 

^^(7! • 72- cos 6) 
47r 2n 

X 47r J ^°& LlLl J rf?7iai^(ii,77i) n(7i,77i,x(?7i,7i)) 

x ]j~ J d\ogL 2 L 2 J dr] 2 a 2 K(L 2 ,ri2) n(L 2 ,7]2, x(r7 2 , 72)) (2-16) 

where a\ = a(rji) and a 2 = 0,(7)2). 

Theoretically we can only predict the ensemble average of fluctuations. To obtain the ensemble aver- 
age of the two- point correlation function, we need the ensemble average (n(Li, 771, Xi)n(L2, ?72, x 2 )) from 
equation ( [2.16] ). In case of point-like sources, it is well known that the ensemble avera ge o f a product of 
the density field has three terms, i.e., homogeneous term, clustering term and shot term |l9| , po| . As we are 
dealing with the density field as a function of luminosity and time, the relation (n(L\, 771, Xi)n(L2, T}2, X2)) 
is not trivial in a strictly sense. However, assuming that the luminosity of sources is statistically inde- 
pendent of their positions relative to the other sources, we set the following relation |l2| , |l3| ]: 

(n(ii,r7i,Xi)n(L 2 ,772,x 2 )) = n(Li,rn)n(L2,ri2) + n(Li, ■q 1 )n(L 2 , f72)£(?7i, t? 2 ; x x - x 2 ) 

+n(L 1 , m )8Qo S L 1 - logL 2 )(5 (3) (xj - x 2 ), (2.17) 



where fi(L,r]) is the spatially averaged quantity of n(L,rj,x). The first term of r.h.s. of (2.17) is the 
isotropic background component, the second term describes the clustering of the sources, the third term 
is called the shot term, which arises from the discreteness of the sources [|l9| . Then the ensemble average 
(C(6)) has three terms, i.e., isotropic background term, clustering term and shot noise term, which we 
write as (C{8)) = (C lso ) + (C ch {8)) + (C SN {8)). 
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A. Isotropic background 



First we consider the isotropic background term. From ( 2.16 ) and the first term of r.h.s. of ( 2.17 ), we 
easily get 



(C IS °) ={hj dlogLL £° dr)aK(L,r,)n(L,rif) = (/«»)' 



(2.18) 



where we used 



^; i ^;^ (7i-72-cosg)=l. (2.19) 



B. Fluctuations from Source Clustering 



The correlation function of the source clustering term is 

(c{e)^) = j f^L^s^.y,- cose) 



X J dlogL 1 L 1 J di] 1 a 1 K(L 1 ,r] 1 ) n(L ll r] 1 ) 



1 



Vc 



X 47r J < ^°^ L ' lL2 J d V2a2K(L 2 ,ri2) n(L 2 ,r]2) 
x £(?ft, 772;xi0?i,7i) -x 2 (?72,72))- (2-20) 
To relate the correlation function of X-ray sources with that of CDM, we write 

£(771,7725x1 -x 2 ) = &x(?7i)&x(r72)£(J7i,?72;xi -x 2 ) C dm, ( 2 -21) 

where bx(v) i s the bias factor^], the CDM density correlation function is £(771, 7725x1 — X2)cdm = 
(5 c (?7i, Xi)5 c (t72, x 2 )\, and <5 c (t7,x) is the CDM density contrast. While the nonlinearity of the CDM 
density contrast may be effective, however, we only consider linear theory in this paper. In this case 
£(771, 7725x1 — X2)cdm is proportional to _Di(t7i)Z?i(?7 2 ), where Di(rf) is the linear growth rate normalized 
to unity at present. 

We next consider to express the CDM density contrast in terms of the density power spectrum. In 
order to express Gaussian random process of the CDM density field, we introduce probability variables 
a-ujim, which satisfy, 

(a-Loim) = 0, (a^imaZ'i' m ') = <K W - u')8u>5 mm >. (2.22) 

For the condition of reality of density field we assume a* Jm = a u i— m . Then we can write the CDM 
density fluctuation fields as 



x) = / rfwVo^^fijjLtatx), (2.23) 
Jo 1 



' The bias factor could be a function of luminosity L in general. However we ignore the dependence of luminosity 
for simplicity. 
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where 5 U (77) is the root-mean-square of the coefficient associated with the orthonormalizcd harmonics in 
the hyperbolic universe 



3W(x) = (-K)V 4 W^i( X )Y lm (n) 



T(iuj + 1 + 1) 



y/sinh X 



which satisfies 



(-X)- 3 / 2 sinh 2 xdxdtoyulmfa n)X>Vm>(x, = S(lu - w')6 u ,6r, 



(2.24) 



(2.25) 



where T(z) is the Gamma function and P£(z) is the Legendre function. Here uj is a non-dimensional 
wave number and K(u> 2 + 1) is the eigen- value of the scalar harmonics 3^>Zm(x) in the hyperbolic universe 
(U) Then we have 



<(5 c (r?i,xi)4(7 ?2 ,x 2 )) = V / 

lm J ° 



duj5 u (r]i)y u i m (xi,0,i)S u (r)2)y* lm (x2,0,2). 



From a straightforward calculation using the relations 

2/ + 1 

<5(7i • 72 ~ cosfl) = ^2 — 9 — P '(7l ' 72)fl( COs6 '): 



Pidi ■ 72) 



47T 



2/ 



m=—l 



where Yi m (fl^) is the spherical harmonics in an unit sphere, equation ( 2.20| ) reduces to ^ 

<C^ CL >=£^Q CL P*(co S 0), 



with 



do;(-if) 3/2 (— / dlogLL / drpK(L,ri)n(L,Ti)bx(v)Su(v)W u i(&Ti) 



(2.26) 

(2.27) 
(2.28) 

(2.29) 
(2.30) 



where Arj = ^f—K{rjQ — rj) and W u i{x) is defined in equation (2.24). This expression is rewritten as Jl5j 



C7, CL = 



dkk 2 M, 



1 

47T 



where M; = (fc 2 - if) • ■ ■ (fc 2 - l 2 K)/(k 2 -K) l ,k = ^/^Kuj, and 

V^ + in^^-^^oshx) 



V smn x 



(2.31) 



(2.32) 



In the limit of flat universe, i.e., K — ► 0, the function X^(x) reduces to the spherical Bessel function 

MM. 



* Throughout this paper we consider an ideal detector with infinite small angular resolution. The effect of a 
finite beam width and a limited observational area can be taken into account by multiplying a window function 
in the similar way as the case of CMB temperature anisotropics. 
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C. Random Fluctuations 



Finally we consider random fluctuations from the shot term. The ensemble average of the two-point 
correlation function is 



(C(0) SN ) 



d f fld % (7i-72-cos( 
47r 2n 



x— / dlogiiLx / d?7iaiJC(Li,77i)— / d\ogL 2 L 2 / dr/ 2 a 2 K(L 2 ,r] 2 ) 



4ir 



4tt 



x n(Li,?7i)(5(Li - L 2 ) £ (3) ( x i("i>7i) - x 2 (772,72)) ■ 
As the delta function is expressed 

^ - X 2 ) = (_ X) 3/2 ^(Xl "X2) (5(2)(fi _ . }) 

sinn xi 



then equation ( 2.33| ) leads to 



c sn w= ^2T^1 c sn P;(cos( 



(2.33) 



(2.34) 



(2.35) 



with 



Cf N - / dlogLL 2 / dr?( —aK(L, v ) I f»(L,o)- 



47T 



sinh A?/ 



(2.36) 



where A77 = \f—K(r]Q — rf). Here we used equations (2.27) and (2.28). This spectrum is independent of 
I, which implies that this fluctuation is a white noise type. And this means (C(#) SN ) oc 5(6). 



III. FLUCTUATIONS IN A SIMPLE MODEL 



In this section we discuss characteristic behavior of fluctuations by considering a simple evolution model 
of X-ray sources. The evolution and population of X-ray sources at high redshift are less well understood. 
Here we consider a model of a single class of sources, i.e., the case that the sources have same luminosity 
profile at a cosmic time. In this case the luminosity function is given by 

n(L,rf) = n(a)5(\ogL-\ogC), (3.1) 

where C is the luminosity of a source and n is the averaged number density of sources. We assume 
that the sources are distributed to the redshift z max . Furthermore we assume power- law evolution of the 
source number density and the luminosity: 

n(a) = n a- d , C = L a~ e , (3.2) 

where no and Lq are constants. We further assume the power-law frequency dependence L v oc v~ a , then 
we have 

K(L,r,) = a(rf) a . (3.3) 

As for the bias, we take into account the epoch-dependent bias by setting pj,|24[ 

bx(v) = bxo + z[b xo ~l}, (3.4) 

where z is the redshift and bxo is a constant. Although this simple formula might be applicable only in 
the Einstein de Sitter universe, we use it even for other cosmological models for simplicity. 
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A. Fluctuations from Source Clustering 



The above simplification leads to the following simple expressions for the isotropic background and the 
angular power spectrum of source clustering: 



where r\i is the conformal time at the redshift z max - Note also that the relative amplitude of fluctuations 
Cj CL /(/(°)) 2 is independent of tlqLq. The evolution of the X-ray sources are less understood. Let us 
consider the case that the volume emissivity evolves with the time, and we set p = d + e = 3 for simplicity 
||. We adopt this value and a = 0.4 in the rest of this paper 

The angular power spectrum is shown in Figure 1 for various cosmological models. In Fig. 1(a) and 
1(b), the three panels correspond to the different cosmological models, a 'standard' CDM model (SCDM), 
a CDM model with the cosmological constant (ACDM) and an open CDM model (OCDM) Q. In each 
panel three lines stand for the different parameters for removal of bright sources. For simplicity we 
take this effect into account by specifying the distance of redshift z m j n instead of specifying S c . We 
assumed that the sources of < z < z m i n were removed. In the figure three lines correspond to the cases 
Zmin = 0, 0.02, 0.1. Fig. 1(a) is the case of no bias by setting bxo = 1 in equation (|3.4|). In Fig. 1(b) 
we chose bxo = 1-6 as a case of the epoch- dependent bias, where the bias factor becomes significant at 
zk,l. The CDM density power spectrum is normalized as erg = 1. 

First notable feature is that the amplitude of fluctuations strongly depends on z m i n for low multipoles, 
i.e., large-angle scales. As z m j n controls distance to which nearby sources are removed, the decrease of 
the low multipoles is due to the removal of the sources. This implies that the nearby sources dominantly 
contribute to the low multipoles Jl5| ]. On the contrary, the high multipoles do not depend on z m j n . 

In order to understand dominant contributors of sources to fluctuations in more detail, we show the 
amplitude of the multipoles I — 1, 10, 30, 100, as a function of z m ; n in Fig. 2. The large value of z m j n 
taken in this figure might be unrealistic from the observational view point for an all sky survey. However 
this figure is instructive to understand what sources are dominant contributors to the fluctuations of 
each angular size. The steep decrease of the low multipole at z m j n ^ 0.1 indicates that low redshift 
sources sensitively contribute to the large- angle fluctuations. On the other hand the higher multipole 
(/ = 100) does not sensitively depend on z m j n and it is not almost affected by the removal of the sources 
for z m j n 2. This means that the high redshift sources of z <> 1 are the dominant contributors to the 
small-scale fluctuations I <> 100. It is shown that the amplitude of the fluctuations increases for z m ; n <^ 2. 
The width of the range in which sources are distributed becomes thin when z m i n becomes near z max . 
This effect makes the amplitude of the fluctuations increase since the fluctuations suffer less significant 
thickness damping, that is contributions from sources at different redshifts cause the cancellation. We 
think that the increase of the fluctuations at z m j n 2 in Fig. 2 are due to this effect. 

Comparing Fig. 1(a) and 1(b), evolution of the bias derives significant difference on small scales. For 
the low multipoles, since the nearby sources of low redshift dominantly contribute to the fluctuations, 
the (low multipole) amplitude scales by the bias factor at present epoch bxo- However the amplitude on 
small scales is significantly increased for the epoch-dependent bias model. This is understood that the 
high redshift sources are the dominant contributors to the small-angle fluctuations and the amplitude is 
increased due to the large value of the bias factor at high redshift. 

It is notable that the dependence on the cosmological models is quite weak even on small scales. Because 
we have employed eg normalization, the shape of the density power spectrum of the OCDM model is 
almost same as the ACDM model. The difference between these two models is small at a fewx 10% level. 
However the amplitude of fluctuations of the OCDM model is slightly larger than the A model on small 
scales. This would be understood as follows. When we see a same physical size in the open universe, the 




(3.6) 



(3.5) 
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angle size becomes smaller compared with the case of the flat universe due to the curvature of background 
geometry. Therefore the large-scale fluctuations in the open model shift to the smaller scale. This is the 
same situation as the case of CMB anisotropics. 



B. Shot Noise v.s. Source Clustering 

In this subsection we discuss the shot noise fluctuations comparing them with the fluctuations due to 
the source clustering. The shot noise fluctuations are random fluctuations originated from discreteness 
of the sources. The spectrum is the white noise type and the angular power spectrum does not depend 



on Z. In our simple model, equation (2.36) reduces to 

SN = r d j a i-e +a y a - d _K 

(47r) 2 J w 'V J sinh 2 Ar7 V ; 

Assuming rj c ~ 770 and (770 — T) c )/r]a C 1, we approximately rewrite it as 



n^Ll r Vc dn ^ npLl 1 
(4tt) 2 J n . (ry - v) 2 ~ (4tt) 2 770 - Vc 



Sc^ , , U , 2 , (3.9) 



where the second equation is derived from (2.15). Then we have 



Cf N ^ S c ^n ^-y 2 . (3.10) 

This approximation is valid for the case that the removed bright sources are located at low redshift z <C 1. 
Using the value = 5.2 x 10" 8 erg s _1 cm _2 str _1 in 2—10 keV band |2(J, we have 



°i N ,J S c \ 1/A f L \ 3/ V n N 1/2 



1(0) L2 ( v erg s-icm-V \3 x 10 43 erg s" 1 ) \3 x 10- 5 Mp C - 3 J ' 

From ( |3.9[ ) and 770 — 77 ~ z/Hq, which is obtained from the Friedmann equation, we estimate the redshift 
z min t° which the bright sources are removed as 

S^) ""(ax 10 43 crg s- 1 ) Gx lO-Hergs-icm- 2 ) ' (3 ' 12) 



As we see from ( 3.10| ), the shot noise fluctuations are controlled by the flux cutoff limit S c . When S c 



is decreased the shot noise fluctuations reduce too. As we saw in the previous subsection, the amplitude 
of the low multipole fluctuations due to the source clustering decreases when S c is decreased. However, 
the higher multipoles are constant even when S c is decreased. The epoch-dependence of the bias factor 
makes the amplitude of small-scale fluctuations significant. Now let us compare these two fluctuations, 
the source clustering and the shot noise. In Fig. 3 the fluctuations due to source clustering are compared 
with the shot noise fluctuations as a function of S c . In this figure we fo cuse d on the multipole of / = 100. 



The upper panel is the case of no bias by setting bxo = 1 in equation (3.4). In the lower panel we chose 
bxo — 1.6. The large difference between these two panels, which is due to the high bias factor at high 
redshift, implies that the bias is a very important factor besides the evolution of the luminosity function 
of sources when predicting the fluctuations due to source clustering. 
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IV. SUMMARY AND DISCUSSION 



In this paper we have formulated equations to investigate the X-ray background fluctuations in various 
cosmological models in a systematic way. We have made full use of the harmonic expansion method in 
terms of multipole components in both open (hyperbolic) and flat universes to obtain the angular power 
spectrum of fluctuations. The fluctuations due to the source clustering and the shot noise can be predicted 
if the evolution of the luminosity function and the bias factor are given. We have calculated the angular 
power spectrum for the X-ray background fluctuations in various cosmological models based on a simple 
X-ray source model. As Lahav, et al. pointed out the fluctuations predicted in the simple model 
strongly depend on the evolution parameter of sources, and rather weakly depend on the cosmological 
models. Moreover the large-angle (low multipole of I & 10) fluctuations are strongly affected by the flux 
cutoff parameter for removal of nearby bright sources |l5|,|| • This is because the nearby sources of low 
redshift are the dominant contributors to the large-angle fluctuations. Therefore when one makes the 
comparison between the theoretical prediction and the observational data on large angular scales, the 
procedure of removing bright sources has to be very carefully taken into account. 

In this paper, we have developed the formulation which is useful to calculate the angular power spec- 
trum. In the observational point of view, however, we may not have to derive Ci in small scales ||||. The 
description by (C(#)) is usually used when one compares the theoretical prediction with observational 
data (e.g., JT^,0). Of course, we can easily derive (C(6>)) from Ct by equation ( 2.29 ) with a window 



function. Moreover, our approach is useful to understand the physical mechanism which determines the 
behaviors of fluctuations on various scales. This formulation is not limited to the X-ray background fluc- 
tuations but it can be easily generalized to calculate the clustering of point-like sources at high redshift. 

In predicting the X-ray background fluctuations the evolution of luminosity function is an important 
factor. Our analysis shows that the evolution of the bias factor is also a very important factor especially 
on small scales. This is because the small-scale fluctuations contain information of clustering in the high 
redshift universe. In order to use the X-ray background fluctuations as a probe of large scale cosmological 
density fluctuations at high redshift, the knowledge of the evolution of the bias factor and the luminosity 
function is necessary factor, because the dependence of cosmology on the X-ray background fluctuations 
is rather weak. Conversely, which suggests that the X-ray background is a good probe for the bias 
mechanism of high redshift sources |2^j2^] when the cosmological parameters are fixed. The future X-ray 
mission projects will give us fine solution for the evolution of X-ray sources. In that case the X-ray 
background will be a possible probe to investigate the clustering and formation process of X-ray sources 
in the high redshift universe. 
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FIGURE CAPTION 



Fig. 1 — The angular power spectrum of fluctuations due to source clustering for various cosmological 
models. The cosmological parameters are taken as ft- = 0.5 and fig = 1 for the SCDM model, and h = 0.7 
and fig = 0.3 for the ACDM and OCDM models. The parameters p = d + e = 3 and z m ax = 3 are taken. 
In each panels three lines correspond to z m ; n = 0, 0.02, 0.1, respectively. Fig. 1(a) is the case of no bias 
by setting bxo = 1 in equation ( |3 . 4[ ) . In Fig. 1(b) bxo = 1.6 is chosen as a case of epoch-dependent bias. 
Note that the area under the lines do not directly describe the amount of observed fluctuations since 



y/£(2£ + l)Ct is the one in case of the logarithmic interval of x-axis (see Eq. 2.29 ). Roughly speaking, 
we need to multiply £. Therefore there are larger fluctuations on smaller scales. On the other hand, the 
contribution from the shot noise does not depend on £ in these figures. 

Fig. 2 — z m i n dependence of multipoles for the ACDM model. The four panels correspond to the multipoles 
I = 1, 10, 30, 100, respectively. The solid line represents the case bxo = 1 and the dashed line does the 
case bxo = 1-6. The model parameters are same as that in Fig. 1(a). 

Fig. 3 — Comparison of fluctuations due to source clustering and shot noise term as a function of flux 



cutoff limit S c . The upper panel is the case of no bias by setting bxo = 1 in equation (3.4). In the lower 
panel bxo — 1-6 is chosen as a case of epoch-dependent bias. Each panel shows the multipole I = 100. 
The solid line shows the clustering fluctuations and the dashed line does the shot noise fluctuations. For 
the clustering fluctuations we used the SCDM model with h — 0.5 and fig = 1. Here we normalized the 
CDM density fluctuations by as = 0.5. In order to calculate the shot noise term, we need to specify the 
values of e and d, separately. Though there are many uncertainties in the evolution model of the sources, 
let us take the values e = 1.4 and d — 1.6. In th is case we have 1 + a — e — 0, which allow us some 
analytic calculations when solving equation ( 2.15| ). We also used the parameters Lg = 3 x 10 43 erg/s, 
n = 10" 5 / M P c3 - We set z max = 3.2 in order to be = 5.2 x 10~ 8 erg/s/cm 2 |§. In this case of 

the parameters the estimation for the shot noise fluctuation is ^Jcf^/I {0) ~ 0.7S* c 1/4 /(erg/s/cm 2 ). The 

dashed line shows the result obtained by calculating Eq. fl3.7| ). The dotted line is the result obtained by 
extrapolating the estimation by Lahav et al. 0]. 
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